A Hermite-Pade" perspective on Gell-Mann— Low renormalization group: an 
application to the correlation function of Lieb-Liniger gas 
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While Pade approximation is a general method for improving convergence of series expansions, 
Gell-Mann-Low renormalization group normally relies on the presence of special symmetries. We 
show that in the single-variable case, the latter becomes an integral Hermite-Pade approximation, 
needing no special symmetries. It is especially useful for interpolating between expansions for small 
values of a variable and a scaling law of known exponent for large values. As an example, we extract 
the scaling-law prefactor for the one-body density matrix of the Lieb-Liniger gas. Using a new result 
for the 4th-order term in the short-distance expansion, we find a remarkable agreement with known 
ab initio numerical results. 
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Quantitative information very often comes in the form 
of a series expansion in a dimensionless variable [l| . Of- 
ten one knows only the first few terms; but even if one 
knows many, the raw series is usually useful only for a 
limited range of values of the expansion parameter. One 
may therefore turn to various methods for improving the 
convergence properties. This is justified if one expects 
that the solution, in the region of interest of the expan- 
sion variable, is an analytic continuation of the solution 
in the region in which one is performing the expansion. 
The best-known method is the Pade approximation Q 
and its various Hermite-Pade generalizations [3J. 

On the other hand, when extrapolations from one pa- 
rameter region to another were called for within particle 
physics, the solution was the Gell-Mann-Low renormal- 
ization group (GML-RG) 0,(1. This method is normally 
not comparable to the Pade approximation because it 
takes as input not only the known terms of the expan- 
sion, but also certain known nontrivial symmetries in the 
problem, an approach that was later formalized and ex- 
tended to the study of differential equations using Lie 
group-theoretic methods @. We should also mention a 
related body of work in Ref . Q • 

We analyze the mechanism of operations of GML-RG 
in the single-variable (SV) case. We show that the SV 
GML-RG requires no special symmetry and is in fact an 
integral Hermite-Pade approximation Q, a connection 
that remains largely unexplored. The method is appli- 
cable to any series expansion, but is especially good for 
interpolation between a power-series expansion (for small 
values of a variable) and a scaling law (for large values of 
the variable) , where the scaling law exponent is a known 
continuous function of a parameter, and one seeks, e.g., 
an approximation for the prefactor. As a case in point, 
we treat the density matrix of the Lieb-Liniger gas l9|] , an 




object of considerable interest to both cold gases 
and integrable systems [l2| communities. 

The GML-RG suggests itself very naturally in prob- 
lems where a physical quantity at one value of a param- 



eter is given as an expansion in terms of the value of the 
same quantity at another value of the parameter, a situ- 
ation usually caused by regularization and renormaliza- 
tion A simple single-parameter example is T(fcb) = 
T{k a ) + 2 f iln(k b /k a )T 2 (k a ) + (2^ 2 ln 2 (k b /k a )T 3 (k a ) + 
• • • , which is the Born series for the scattering prob- 
lem for the 2D delta- function potential T(k) is the 
complex-valued function one is trying to determine (T- 
matrix), \i > is a constant (mass), and the expansion 
holds for all k a , k b > (momenta). The expansion would 
not seem to be useful unless kb/k a ~ 1; however, we may 
try a "divide and conquer" approach: instead of going 
from k a to k b in a single large step, we use many smaller 
steps, k a = kg — > fci — > ki — > ■ • • — > fc/v = k b , for each 
of which kj+i/kj ~ 1. If k n+ i/k n — 1 + e, we may let 
e — » 0, N — > oo so that (1 + e) = const. = k b /k a . 
We obtain an equation of the form kT'(k) — f3(T(k)), 
the Callan-Symanzik equation (CSE) for the problem, 
whose solution is the RG-improved solution of the prob- 
lem. In this example, (3(T) = 2/j.T 2 at all Born or- 
ders. Thus, starting with the first nontrivial Born or- 
der, the RG treatment here produces the exact result, 
l/T(k b ) = 1/T(k a ) - 2fj,]n{k b /k a ). Now note that the 
CSE may be obtained by applying d/dk b \ kb ka ^ k to 
both sides of the expansion. Thus, to get the RG- 
improved solution, we differentiate the original expan- 
sion and then promptly integrate it, which raises the 
question: why do we not end up right back where we 
started? More formally, let T(k) be the quantity of 
interest; then for some r, r = k r T is dimensionless. 
The exact solution of the problem may be written as 
i~(k b ) = H(r(k a ), k a , k b ) for some function H. Dimen- 
sional analysis then says that there must exist a function 
G such that r{k b ) = G [r(k a ), k b /k a ] (scale invariance or 
self- similarity of r). Applying d/dk b \ kb fc a ^, we get the 
CSE, kr'{k) = /3(r(fc)), where (3{x) = G^ ' x \x, 1). This 
is easily integrated, |^ = exp dt, and we see 

that if we start with the exact dependence of r(k b ) on 
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r(fc Q ) and kb/k a (i.e. exact G), the RG procedure must 
give the same exact dependence back. On the other hand, 
if we start with a truncated series expansion of G (a finite- 
degree polynomial in r(fc a )), then the RG does not give 
the same finite-degree polynomial back. But how does 
the RG "know" that the polynomial was not itself the 
exact G? The answer is that not every function of two 
variables can play the role of exact G, because an exact 
G must satisfy a special group property ||. Namely, if 
we wish to relate r(kb) and r(k a ), we can do it directly, 
r(kb) = G(r(k a ), kb/k a ), or we can go through an inter- 
mediary, r(k a ) — > r(fc c ) and then r(fe c ) — ► r(fc(,). The re- 
sults will agree only if G satisfies G(x, ab) = G(G(x, a),b) 
for all a, b, and x. This is a group property, because it 
allows us to introduce a (Lie) group with elements Q a by 
defining group action on numbers through Q a x — G(x,a). 
Group multiplication then corresponds to composition, 
GbGaX — G(G{x, a),b) = G(x, ab) = Q a bX- The generator 
gives the /3-function: Bx = Qa\ a=1 x = (3(x). 

The analysis that follows is straightforward, though it 
seems that it has not yet been reported in the literature. 
First, if G(x, z) is to satisfy the group property, it can- 
not be a polynomial of finite degree in x, because then 
G(G(x,a),b) would be a polynomial whose order is the 
square of that of G{x, ba) . Second, for any function of two 
variables F(x, z) (regardless of whether is satisfies the 
group property or not) we may define /%) = F^ ' ^(x, 1). 
The solution of the CSE is a function F(x, z) which sat- 
isfies the group property, has the same /3-function as 
F(x, z), and is equal to F(x, z) if and only if F(x, z) itself 
satisfies the group property. Next, suppose F(x, z) is a 
truncated expansion in x of some true G{x, z), so that it 
does not satisfy the group property; the key question is 
why t(x) — F(tq, x/xq) is so frequently a better approxi- 
mation of Tt rU e(x) = G(tq , xj xq) than t(x) = F(tq,x/xq) 
is. At first glance, the answer may seem obvious: be- 
cause Tt rue (a;) is expressed in a self-similar form (through 
G), its approximation should also be represented in a 
self-similar form; t(x) is so represented (through F), 
while t(x) is not. However, this cannot be right, be- 
cause all functions of single variable are representable in 
a self-similar form, at least on intervals where they are 
monotonic: f(x 2 ) = /(ffxl) = f( S f-\ f{xi) )). (Prov- 
ably, a function G satisfies the group property if and 
only if there is a single-variable function / such that 
G(x,z) = f(zf^ 1 (x)).) When this procedure is applied 
to t(x) — F(tq,x/xo), we may call it "an upgrade to a 
self-similar form through an exact inversion of the ap- 
proximate expression." However, instead of exactly in- 
verting the polynomial that is F(tq,x/xo), we may also 
think of it as a truncation of an infinite series and per- 
form an order-by-order series reversion; then we com- 
pose the truncated series and its reversion according to 
f{zf~ 1 {x)). This way, any expansion in a single vari- 
able can be represented as a truncated self-similar ex- 
pansion, which is then treatable by the GML RG (see 



below for a simpler method). However, the question of 
why t(x) is "better" than t(x) remains. To proceed, 
let G be a function satisfying the group property; let 
z = X2/X1 so that t{x2) — G(t(xi), z) and assume G 
has an expansion in powers of t\ = t(x\): G{t\, z) — 

fi(z) n + h( z ) T i + M z ) T i H ■ Imposing the group 

property and equating terms of like orders in t\, we 
obtain a recursive sequence of functional equations: 
h(ba) = h(b)h{a); Mba) = A 2 (6)/ 2 (a) + f 2 (b)h(a); 
f 3 (ba) = f!(b)h(a) + 2/ 1 (6)/ a (6)/ a (o) + h(b)h(a); . . . 
We get fi{z) = z r , where r is a real number, and 
f 2 (z) = a 2 / 2 .2, h{z) = a\ / 3 . 2 + a 3 / 3)3 , etc., where 
ha = z -{z r ~ 1), h,2 = f(z r - l) 2 , etc. (If r = 0, 
we have / 2 .2 = kiz; / 3 . 2 = In z; / 3 3 = lnz, etc.) Here 
a-2 = f% (1), a 3 = / 3 (1), etc., are arbitrary constants. It 
will be convenient to introduce an expansion in "rows:" 

= z r Tl +0+0 +0 + .» 

+ aih,1 T \ + a 2f3,2 T l + »2/l,2 T l + ■■■ 

+ 03/3,31"]* + a 3 al / 4>3 rf + ■■■ 

+ tt4 /4,4 Tj 1 + 

Note that f n { z ) T i is the sum of the entries of the nth 
column; the nth order of perturbation theory gives the 
first n columns. Suppose we knew the expansion to the 
second order. We see that if / 2 (z) is nonzero, then a 2 is 
nonzero, and so even though we do not know / 3 (^), we 
can nevertheless conclude that it must at least have the 
term proportional to a\. In fact the entire second row in 
the table must be nonzero. (In the 2D scattering example 
above, all the other rows were zero.) More generally, from 
the first n columns we can deduce the full first n rows. 
What RG does, in effect, is fill out all the rows whose 
entries can be deduced from the known columns, sets the 
other rows to 0, sums the result to infinite order, and 
analytically continues beyond the radius of convergence 
of the resulting infinite series. To see this, note that the 
/3-function is determined solely by the leading terms of 
each row: /?(n ) = r t\ + a 2 rf + a 3 rf + ■ ■ ■ . The solution 
of the CSE with the /3-function truncated to the nth or- 
der then is precisely an Fir, z) whose "row expansion" 
has just the first n rows. In short, the RG treatment re- 
sults in a minimal upgrade to a self-similar form: it only 
includes the rows we know how to reconstruct from the 
perturbation expansion. In contrast, an upgrade through 
an exact inversion of the approximate expression will in 
addition produce spurious rows, artifacts that have noth- 
ing to do with the underlying problem. Note that when 
we use an ordinary perturbation expansion to compute 
actual numerical values, we set to zero all orders we did 
not actually compute; in RG we simply follow the same 
rule, though this time we apply it to the rows. 

There is also a completely different way to think about 
SV GML-RG, which makes it clear that it is a species of 
Hermite-Pade approximation. From the representation 
G(x,z) = f[zf~ 1 (x)] it is easy to show that (3{f}(x) = 
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/ 1 (x)/[f 1 ]'(x). Now consider integrals of the form 
exp f^ b l/g{t) dt. If g(t) — (3{f}(x), the integral is iden- 
tically equal to f^ 1 (xb)/f~ 1 (x a ), and represents the inte- 
grated CSE. We are interested, however, in the situation 
when we only know a truncated expansion for /, T{/}. 
By formal series manipulation, we can compute two ob- 
jects that are asymptotically (for small x) equal to /?{/}: 
(a) /?{T{/}}, and (b) T{f3{f}}. If g is set to (a), namely 
to Tlf-^/Tilf- 1 ]'}, the integral exp f** l/g(t)dt will 
return T{f- X }(x b ) /T{f- 1 }(x a )- Setting x a = and per- 
forming a series reversion then brings us right back to 
where we started, to as it must because every step 

was reversible. However, if g is set to (b), namely to 
T{f~ 1 /[f~ 1 ]'}, that is an irreversible step, and the result 
is something new: an RG approximant. More formally: 
suppose / has an expansion (for simplicity about zero) 
cox a ° + c\x ai + ■ ■ ■ with all Cj ^ and Reao < Reai < 
• • • . Let T n {f}(x) be the truncation of that expansion 
up to and including the x Q ™ term. Now formal series 
manipulations allow one to compute truncated expan- 
sions of various transformations of / (such as / 2 , 1//, 
/', etc.) order-by-order, up to some maximal order 
that is still completely determined by T„{/}(i). Let us 
denote these derived expansions by, e.g., T[/ 2 | T„{/}](x) 
(for example, T[l/sin |T 2 {sin}](x) = ~ + \ x). Now we 
introduce the (order-a„, centered-at-O) renormalization 
group on-the-inverse approximant (RGIA) through 

RGIA„{/}(x) = 

f x 1 

lim T n {f}(e) exp / - T —, — : — ds . (1) 

«-o nlJn 1 P J e T[///'|t„{/}]( s ) V j 

This expression is the key result of the paper. Note that: 
(1) the usual RG result is given by RGIA„{/- 1 }(x); (2) 
usually T[f/f | t„{/}](s) is a polynomial, in which case 
the integral can be evaluated analytically (in terms of 
logarithms and arctangents) once the polynomial is de- 
composed into partial fractions. This requires factoring 
the polynomial, in other words, finding its zeros, which 
is normally the only numerical step in the computation. 
Replacing the denominator in Eq. {T]) by its Pade ap- 
proximation tends to improve results further. 

There is at least one class of functions / for which it 
is obvious that RGIA„{/} must be a better approximant 
for / than T„{/}: the functions whose only singulari- 
ties are poles (of any order) but whose derivatives have 
no zeros. The reason is that for such functions, ///' is 
analytic everywhere. An example is f(x) — tana;, with 
///' = |sin2x. Near the pole at x = n/2, the Taylor 
expansion for / fails as an approximation for /. However, 
the Taylor expansion for ///' remains a good approxima- 
tion for ///', and so RGIA„{/} remains a good approx- 
imation for /. Order-by-order, in this case, at least, it 
also performs better than the Pade approximation. One 
should emphasize that experience thus far shows that the 
RGIA is useful for a much broader class of functions, de- 



spite the fact that there is no clear reason why it should 
be so; but in the case of Pade approximation, it is also 
still not completely understood why it performs as well 
as it does, even after more than a century of research 0]. 

It is however clear that RGIA is a special case of 
the so-called integral Hermite-Pade approximant Q, de- 
fined as follows: one finds the polynomials P(x), Q(x), 
and i?(x), of orders p, q, and r, that satisfy the condi- 
tion P(x)f(x) + Q(x)f (x) + R(x) = 0{ X v +q+r+2 ), with 
P(0) = 1; the more terms one has in the expansion of /, 
the higher the order of agreement, p + q + r + 2, that one 
can demand. Then one sets the right hand side to zero; 
the solution of the resulting equation is the integral ap- 
proximant. We see that RGIA corresponds to R(x) = 
and Q(x) = -/'(0)//(0). Thus, the SV GML-RG is 
an integral Hermite-Pade approximation plus a series re- 
version 0: P^/" 1 ]'^) - l -^p^f- 1 (x) = 0(xP+ 2 ). 
We note in passing that while in the Pade approxima- 
tion community it is believed that best approximants 
should be invariant under the homographic transforma- 
tion x = Aw/(1 + Bw) [l5j |. this invariance demands 
p = q + 2 (and q = r), so that the RGIA is invariant un- 
der this transformation only if p = 2. On the other hand, 
high-energy physics community has accumulated a huge 
amount of knowledge and experience about GML-RG. In 
the future we plan to bring the two research traditions to- 
gether around the point of contact, the RGIA, and begin 
to work out what mutual implications they may have. 

Turning to practical applications of the RGIA, we 
will consider problems where, in addition to the ex- 
pansion of f(x) for small values of x, one also knows 
the exponent of the scaling law for large values of x: 
f(x) ~ C/x a with a known but C unknown. The 
idea is to apply the RGIA to f(x) — x a f(x) (one can 
treat the case where / ~ C(xq — x) 13 at a finite xq in 
a similar way). Problems like this can be approached 
using Pade approximation as well [lil Il6j |. but for 
exponents which are not simple fractions, the Pade 
method becomes inconvenient. For example, to treat an 
irrational exponent, one would have to approximate it 
by a ratio of integers p/q, and then find the approximant 
for [f(x)] q . This becomes particularly unnatural if the 
exponent is a function of a continuous parameter. In this 
case RGIA is clearly a better choice, because even for 
irrational exponents a, p = f / f is a rational function; 
the RGIA will then itself be a continuous function of 
the exponent. As a nontrivial application, we con- 
sider the one-body-reduced density matrix p\(z\; z[) = 
J Q ...J Q dz 2 ... dz N ^(zi, z 2 , . . . , z N )^*{z' l7 z 2 , . . . , z N ) 
of the ground-state Lieb-Liniger gas (a ID gas of bosons 
interacting via interparticle potentia l q -\pS(z)) with 
periodic boundary conditions 0, 0, [HI, 0|. In the 
thermodynamic limit, where the number of particles N 
and the system length L go to infinity but the ID density 
n = N/L is kept constant, L/?i(0; , z) = g\{z)jn is finite 
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FIG. 1: One-body-reduced density matrix of the Lieb-Liniger 
gas for a variety of interaction strengths 7. Thick lines: re- 
sults of an ab initio quantum diffusion Monte Carlo (QDMC) 
simulation; thin lines: an RGIA interpolation between the 
fourth-order short-distance expansion and the long-distance 
scaling law. Not labeled are the nearly horizontal curves for 
7 = 0.001. To the right of the plot is a table comparing 
the values for the prefactor of the long-distance scaling law 
obtained from the RGIA interpolation, and from the QDMC. 



obvious direction would be to see if invariance under no- 
mographic transformation plays any role in GML-RG. 
Finally, we singled out a class of interpolation problems 
for which GML-RG is particularly well-suited; we demon- 
strated the effectiveness of the method in the case of the 
density matrix of the Lieb-Liniger gas, for which we have 
also announced the exact form of the fourth-order term 
in the short-distance expansion. 

We thank R. Shakeshaft, M. Jaffrey, and M. Rigol for 
fruitful discussions, and G. E. Astrakharchik for provid- 
ing us the data from his QDMC calculations. This work 
was supported by the Office of Naval Research grant 
No. N00014-06-1-0455 and the National Science Foun- 
dation grant No. PHY-0754942. 



and its short-distance expansion, nz — + 0, is given as 

gx{z)/n = 1+02(7) (nz) 2 + C3{l) \nz\ 3 + c 4 (j) (nzf-\ . 

Here 7 = [m/fi 2 )giu/n is the dimensionless strength 
of the interaction. In [l(J we showed that €1(7) = 0, 
c 2 (7) = -5 [e 2 (7) -7e 2 (7)] ; and c 3 {j) = ^tMM, 
where e m (^) — [7/A(7)] m+1 /_ 1 x rn g (x\j) dx is the nor- 
malized rnth moment of the "density of rapidities" func- 
tion g(x\~f), which, together with the function A(7), is 
determined through the Lieb-Liniger system of equations: 
2tt 5 (x\j) - [2A( 7 )]/[A 2 ( 7 ) + (y - x) 2 ] g (y\i) dy = 1 
and 7 J j g (xl'y) dx — A(7). We also use a new 
result, 04(7) = 7^/12 - 3e 4 /8 + (2 7 2 + 7 3 ) e' 2 /24 - 
7 62/6 — 7 €2 e;>/4 + 3 e?,/4 , whose derivation will appear 
elsewhere [lTj- For nz — > 00, from Luttinger liquid 
theory [ll| it is known that g\{z)jn = C(j)/ (nz) a ^ 
with 87r 2 a 2 = 662(7) — ^l e 2 + 7 2e 2i C(7) is not known, 
and we extract it using GML-RG. In Figure [IJ we 
compare the results of an RGIA interpolation between 
the short-distance expansion and the long-distance 
scaling law for g\{z)jn with the results of an ab initio 
quantum diffusion Monte Carlo simulation llj, at a 
wide range of interaction strengths 7. In the table to 
the right of the plot, we compare the results for the 
prefactor C(7), finding remarkable agreement. (We 
should mention here the excellent analytical estimate 
of (7(7) from an inspired, though as-yet unjustified, 
extrapolation of Popov theory in Ref. [191].) 

We have analyzed the mechanism of operation of 
single-variable Gell-Mann-Low renormalization group 
(SV GML-RG). We found that it does not rely on any 
special symmetries and that it is in fact is a combination 
of series reversion and an integral Hermitc-Pade approx- 
imation. To the best of our knowledge, such an explicit 
connection between the GML-RG and Pade methods had 
not yet been pointed out, and future work may use this 
connection to explore the implications that the GML-RG 
and Pade research traditions have for each other. One 
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